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Abstract 

The Neogene eolian deposits in tlie Cliinese Loess Plateau (CLP) are one of tine nnost useful continental deposits for 
understanding climatic changes. To decipher Late Neogene paleoclinnatic changes in the CLP, we present a terrestrial 
mollusk record spanning the time interval between 7.1 and 3.5 IVla from the western CLP. The results indicate four stages of 
paleoclimatic evolution: From 7.1 to 6.2 Ma, cold and dry climatic conditions prevailed as evidenced by high values of the 
total number of cold-aridiphilous (CA) mollusk species and by low values of all of the thermo-humidiphilous (TH) mollusk 
indices. From 6.2 to 5.4 Ma, the climate remained cold and dry but was not quite as dry as during the preceding phase, as 
indicated by the dominance of CA mollusks and more TH species and individuals. From 5.4 to 4.4 Ma, a warm and moist 
climate prevailed, as indicated by high values of the TH species and individuals and by the sparsity of CA species and 
individuals. From 4.4 to 3.5 Ma, all of the CA indices increased significantly and maintained high values; all of the TH indices 
exhibit high values from 4.4 to 4.0 Ma, an abrupt decrease from 4.0 Ma and a further increase from 3.7 Ma. The CA species of 
Cathaica pulveraticula, Cathaica schensiensis, and Pupopsis retrodens are only identified in this stage, indicating that the CA 
species were diversified and that the climate was becoming drier. Moreover, the CA mollusk group exhibits considerable 
diversity from 7.1 to 5.4 Ma when a cold, dry climate prevailed; whereas the diversity of the TH group was high during the 
relatively warm, wet interval from 5.4 to 4.4 Ma. This indicates that variations in the diversity of the CA and TH mollusk 
groups were closely related to climatic changes during the Late Miocene to Pliocene. 
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Introduction 

As evidenced by the marine benthic foraminiferal S'^'O record, 
Earth's climate underwent a gradual cooling trend during the Late 
Miocene and Pliocene. This interval witnessed the progressive 
cooling of oceanic deep water, the expansion of permanent ice 
sheets in Antarctica, the occurrence of ice rafted detritus in 
Northern Hemisphere (mostly in north Atlantic), and both 
hemispheres covered by ice sheets after middle Phocene [I]. 
Coincident with ice development in both Polar regions are 
significant ecological, climatic, and tectonic events elsewhere and 
especially around Asia [1-12], and which demonstrate that the 
Late Miocene and Pliocene was an important and complex 
interval that needs to be better understood. 

In East Asia, one of the most important climatic changes is the 
evolution of the East Asian (EA) monsoon. Numerous sedimen- 
tological, geochemical, and paleontological studies, including of 
fossil mammals, mainly from the Chinese Loess Plateau (CLP) and 
the South China Sea (SCS), have contributed significantly to our 
understanding of the EA monsoon changes during the Late 



Miocene to Pliocene. However, the results are inconsistent and 
some are even in conflict. In the CLP, sedimentological evidence 
of changes in sediment grain size and sedimentation rate from the 
Xifeng, Lingtai, Lantian, Xunyi, Luochuan, Jiaxian, Baode, and 
Jingle (Figure 1) Red Clay deposits indicate the occurrence of a 
strong winter monsoon and pronounced aridity in the Asian 
interior during the Late Miocene [8,13-16]. However, the 
geochemical record of Sr/ Sr ratios from the Lingtai Red Clay 
sequence implies a weak EA winter monsoon from 7 to 2.5 Ma 
[17]; and in addition, this interval can be sub-divided into a large 
amplitude and high frequency stage from 7^.2 Ma and a weak, 
stable stage from 4.2-2.6 Ma, as evidenced by records of Zr/Rb 
ratio and mean quartz grain size [18]. Furthermore, grain-size 
records from fluvial deposits in the Linxia Basin (Figure 1) indicate 
that the EA winter monsoon intensified at 7.4 Ma and 5.3 Ma 
[19], the latter datum being significandy different to evidence from 
the Red Clay record. 

In the case of the EA summer monsoon, the magnetic 
susceptibility record from the Xifeng Red Clay sequence indicates 
a weakened summer monsoon during the Late Miocene (6-5.4 
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Figure 1. Location of the Chinese Loess Plateau (CLP), the studied Dongwan loess-paleosol sequence and other Neogene sections 
mentioned in the text. The Dongwan sequence is indicated by a red circle and other sections by white circles. Main cities are indicated by black 
squares. 
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Ma) and a strong summer monsoon during the early Pliocene [8] . 
Additional information about the EA summer monsoon is 
provided by records of pedostratigraphy and iron geochemistry 
from the Lingtai Red Clay deposits. These records indicate that 
the EA summer monsoon was relatively weak from 7.05 to 6.2 Ma, 
strengthened from 6.2 to 5.5 Ma, was very strong from 5.5 to 3.85 
Ma corresponding to the interval of strongest soU development, 
and weakened significantly over the interval 3.85 to 3.15 Ma 
[20,21]. 

However, studies of mammalian fossils from the CLP have 
yielded very different results: Hypsodonty analysis indicates that 
northern China became more humid at 7-8 Ma, coincident with 
the onset of Red Clay deposition in the eastern CLP, and which 
was interpreted by the authors as representing the onset or 
intensification of summer monsoonal precipitation [22,23]. A 
general rise in 8 O values of soil carbonate from the Lantian 
fluvial and Red Clay deposits reflects increased summer precip- 
itation related to the onset and/or intensification of the EA 
monsoon during the Late Miocene to Pliocene [24,25], and which 
is supported by a mammalian faunal turnover event implying a 
marked change to more humid and closed habitats [26]. Recent 
isotopic evidence from fossil mammals and soil carbonates 
indicates a strengthened EA summer monsoon from 7-4 Ma 
[27]. However, the S C values of fossil enamel from a diverse 
group of herbivores and of paleosol carbonate and organic matter 



from the Linxia Basin indicate that C4 grasses were either absent 
or insignificant in the Linxia Basin prior to 2—3 Ma, suggesting 
that the EA summer monsoon was probably not strong enough to 
affect this part of China throughout much of the Neogene [28]. 
However, this result conflicts with the interpretations of Fortelius 
et al. (2002) [22], Kaakinen et al. (2006, 2013) [24,25], Liu et al. 
(2009) [23] and Passey et al. (2009) [27]. 

In the SCS, the EA winter monsoon developed progressively 
from 7 Ma as shown by an increasing trend in black carbon 
concentration and accumulation rate [29]. In addition the EA 
summer monsoon weakened after 7.5 Ma, as evidenced by 
combined planktonic foraminiferal Mg/ Ca and 5 O records [30] , 
consistent with the Red Clay sequences in the CLP. The 
suggestion of a weakened EA summer monsoon is supported by 
geochemical records which indicate that chemical weathering 
intensity decreased from the early Miocene with a rapid decrease 
centered at 7.2 Ma [31]. This result is generally consistent with 
records of radiolarian species numbers and individuals, and 
diversity, from the southern SCS and which suggest a major 
decrease in summer monsoon intensity after 7.70 Ma [32]. In 
contrast, records of clay/feldspar ratio, kaolinite/ chlorite ratio and 
biogenic opal MAR from the SCS suggest that the summer 
monsoon was strong from 7.1-6.2 Ma and remained relatively 
stable from 6.2-3.5 Ma [33]; however, the authors emphasize that 
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their study is a schematic reconstruction which only outlines the 
principal stages but neglects the details. 

A semi-quantitative reconstruction of the Neogene vegetation in 
China indicates that Miocene aridification associated with 
strengthening of the EA winter monsoon is consistent with 
Neogene global cooling, and that the EA summer monsoon did 
not weaken during the Pliocene [34], and this result agrees with 
most of the geological records from the CLP and SCS. However 
another quantitative reconstruction from plant fossil records 
yielded contrasting results, indicating that records of both 
temperature and precipitation from north China exhibit no 
significant difference between the western and eastern regions 
during the Miocene, suggesting that the monsoon climate did not 
commence or intensify at that time [35]. Regional climate model 
experiments also reveal that during the Late Miocene, from 1 1-7 
Ma, the monsoonal climate may not have been fully established in 
various Asian regions, including northern China [36]; and this 
finding is contrary to that of numerous previous studies of the EA 
monsoon which suggest that it was initiated around the time of the 
Oligocene/Miocene boundary [37-39]. 

The foregoing review demonstrates that more work is needed to 
better understand the evolution of the EA monsoon during the 
Late Miocene and Pliocene. In particular, higher resolution studies 
using more sensitive monsoon proxies from key monsoon- 
dominated regions may be one of the best solutions for resolving 
the various inconsistencies and even conflicts regarding the process 
of monsoon development. 

As mentioned above, the CLP (Figure' 1), located to the 
northeast of the Tibetan Plateau, is a key continental region for the 
study of the EA monsoon. Deposition of eolian sediments 
commenced in the western CLP from 22 Ma, as observed in the 
QA-I, QA-II, and QA-III (Figure 1) Miocene loess-paleosol 
sequences which have a basal age of about 22 Ma [37,40], while 
the upper boundary is dated at about 3.5 Ma in the Dongwan late 
Miocene-Pliocene loess-paleosol sequence in the western CLP 
[41]. In contrast, in the eastern CLP, the age of the lower 
boundary of eohan deposits, the Red Clay sequences, is about 7-8 
Ma as evidenced in the Lingtai, Xifeng, Liantian, and Baode 
sequences [15,20,42,43]. However, results from the recently 
reported Shilou Red Clay serjuence indicate that eolian sediments 
were deposited from 11 Ma in the eastern CLP [44]. Despite the 
occurrence of a totally different lower boundary age between the 
east and west CLP, these deposits have great potential for 
deciphering in detail the processes of ecological and climatic 
evolution in the CLP, and by extension in East Asia, during the 
Miocene and Pliocene. As reviewed above, sedimentological and 
geochemical studies have so far contributed a great deal to our 
understanding of climatic changes during the interval of interest 
[6,8,9,13-18,20,21,24,27,37,40^5]; however, to date tiiere has 
been only a limited application of a biological approach to 
analyzing these sequences [22,26,46,47]. 

Land snails are generally the most common and abundant 
fossils in Quaternary loess sequences, and this fact, together with 
the limited degree of success of most of the other paleontological 
studies of Quaternary loess, makes them especially important 
paleoenvionmental indicators for loess deposits. Their occurrence 
in Quaternary loess was first documented in Europe in the early 
1820s [48], and since then they ha\'e contributed significantiy to 
understanding the origin and paleoclimatic evolution of Quater- 
nary loess deposits in Eurasia and especially in the CLP [46,49- 
57]. However, terrestrial moUusks preserved in the Neogene 
sequences in the CLP have not been investigated in detail until 
recendy. Although they have provided crucial paleontological 
evidence for the wind-blown origin of the Neogene loess sequences 



in the western CLP [58,59], the ecological and climatic 
information that terrestrial mollusks may provide has not been 
well deciphered, with the exception of the record from the Xifeng 
Red Clay sequence spanning the interval from 6.2 to 2.4 Ma in 
the eastern CLP [46] . Furthermore, it is unknown how terrestrial 
moUusk diversity varied in the CLP against the background of 
climate changes during the Neogene, since the necessary studies 
have not been performed. In this study, we studied the record of 
terrestrial mollusks preserved in the Late Miocene to Pliocene 
Dongwan section in order to investigate the evolution of 
paleoclimate and terrestrial moUusk diversity in the western CLP 
during the Late Miocene to Pliocene. 

Materials and Methods 

The Dongwan loess-paleosol sequence (105°47'E, 34°58'N) 
[41] is located in the northeastern part of Qin'an County in the 
western CLP (Figure 1). The current climate in Qin'an County is 
mainly controlled by the EA monsoon, with mean annual 
precipitation of about 400-500 mm and mean annual tempera- 
ture of ~10.4°C. Mean July and January temperatures are 
~22.7°C and — 3.4°C, respectively. Vegetation in the region 
corresponds to a semi-arid temperate steppe [60]. 

The Dongwan sequence, located in the western CLP, is the first 
counterpart of the Red Clay sequences in the eastern CLP [41]. 
The section is about 73.7 m thick, and is composed of 84 
distinguishable loess-paleosol couplets. The chronology of the 
sequence has been established by Hao and Quo (2004) [41] using 
magnetostratigraphic and micromammalian studies (Figure 2). 
First, Hao and Quo (2004) [41] ascribed the approximate age of 
the sequence to the Late Miocene to Pliocene using micromam- 
malian fossil teeth sampled from 20 depths of the Dongwan 
sequence. Second, they defined a series of magnetozones based on 
319 oriented samples collected at 20 or 25 cm intervals and 
correlated them with the Geomagnetic Polarity Timescale [61]. 
Finally, they established a chronology using paleomagnetic 
reversals as age controls and interpolated in between them using 
the magnetic susceptibility age model developed by Kukla et al. 
(1990) [62]. Their results yield an age duration of 7.1 to 3.5 Ma 
for the Dongwan loess-palaeosol serjuence [41]. 

A total of 310 moUusk assemblages were collected from the 
Dongwan sequence using continuous 20-cm-thick samples; howev- 
er, some intervals were sampled at intervals of 10-50 cm, according 
to lifhological changes [58] . About 30 kg of sediment were obtained 
for each sample. In the field, we progressively broke each sediment 
sample into smaller pieces of about 0.5 mm in diameter, at the same 
time collecting all available shells and visible broken pieces. No 
necessary permits for the described field investigations were needed. 
In the laboratory, we attempted to restore any broken shells, and 
then identified and counted them under a binocular microscope. All 
of the identifiable moUusk remains were considered in the individual 
totals using the method of Puissegur (1976) [63]. All of the moUusk 
sheUs are stored in the Institute of Geology and Geophysics, Chinese 
Academy of Sciences. 

For each moUusk assemblage, the numbers of species (S) and 
individuals were counted and a diversity index was calculated for 
all species, thermo-humidiphUous (TH) species, and cold-aridiphi- 
lous (CA) species in order to investigate changes in terrestrial 
moUusk populations and in different ecological groups. We used 
the most widely applied Shannon index [64], sometimes referred 
to as the 'Shannon-Weaver' index and sometimes as the 
'Shannon-Wiener' index by researchers, in order to calculate 
the values of diversity, H (S), of the total, TH, and CA species, as 
foUows: 
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Figure 2. Chronology of the Dongwan loess-paleosol sequence (modified after Hao and Guo, 2004) [41]. 
doi:1 0.1 371/journal.pone.0095754.g002 




where n, is the density measure (in this case the mollusk 
individuals) of the i-ih species (/ varying between 1 and n); S is 
the number of species in the sample, and 



!=1 

The theoretical maximum (i/max) of diversity in any sample is 

expressed as 

Hmax=i0g2S 

EquitabiUty (or evenness, E) is expressed as 



^max 

These indices have been applied to European and North 
American Quaternary terrestrial mollusk assemblages [65-67]. 



Results 

MoUusk fossils are relatively abundant in the Dongwan 
sequence with significant concentrations at ~2 m, 10 m, 30 m, 
50 m, and 70 m depth. Amongst the total of 310 samples, 298 
yielded 16439 mollusk individuals and 12 samples were barren. 
The maximum count reached 1121/30 kg at 54 m depth 
(Figure 3). Variations in total moUusk individuals parallel 
fluctuations in magnetic susceptibility [58], indicating that 
pedogenic processes such as carbonate dissolution did not affect 
the preserved assemblages [58]. 

A total of 24 mollusk species were identified in the Dongwan 
section. They are all terrestrial taxa and consist of Gastrocopta 
armigerella (Reinhardt, 1877), Gastrocopta sp., Punctum orphana 
(Heude, 1882), Punctum sp., Metodontia huaiensis (Crosse, 1882), 
Metodontia yantaiensis (Crosse et Debeaux, 1863), Metodontia 
heresowskii (MoeUendorff, 1899), Metodontia cf. huaiensis, Metodontia 
cf. yantaiensis, Metodontia cf heresowskii, Metodontia sp., Kalklla sp., 
Macrochkmys sp., Opeas sp., Cathaica sp., Cathaica pulveratrix (Martens, 
1882), Cathaica pulveraticula (Martens, 1882), Cathaica schensiensis 
(Hilber, 1883), Cathaica placenta (Ping et Yen, 1933), Papilla aeoli 
(Hilber, 1883), Pupilla grahaui (Ping, 1929), Pupilla sp., Vallonia sp., 
and Pupopsis retrodens i'Martcns, 1879). 

All of these species, except Pupilla grahaui and Pupopsis retrodens, 
have been previously identified in the Chinese Quaternary loess- 
paleosol sequences and most of them have modern representatives. 
For example, Cathaica pulveratrix, C. pulveraticula, C. schensiensis, and 
Pupilla aeoli are the most common species that prefer living in 
relatively cold, dry environments, and are presently distributed in 
northwestern China. They have been regarded as indicative of a 
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Figure 3. Variations in terrestrial moliusics versus deptKi in the Dongwan sequence. Lithology and total mollusk individuals are from Li et 
al., (2006a) [58] and Li et al. (2008) [47]. S(CA)-Total number of species of the cold-aridiphilous (CA) mollusk group. H(CA)-Diversity of the CA mollusk 
group. S(TH)-Total number of species of the thermo-humidiphilous (TH) mollusk group. H(TH)-Diversity of the TH mollusk group. S(total)-Total 
number of species of the total mollusk group. H(total)-Diversity of the total mollusk group. E(total)-Equitability of the total mollusk group. 
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strengthened winter monsoon [49-58]. Conversely, Metodontia 
hmiensis, M. yantaiensis, M. heresowskii, Gastrocopta armigerella, and 
Punctum orphana are species living in a warmer and more humid 
environment, and are distributed mainly in the southeastern part 
of the CLP, where the warm, moist summer monsoon brings 
sufficient precipitation [49-58] . Therefore based on their present 
requirements of moisture and temperature, as well as their modern 
geographical distribution, the Dongwan terrestrial moUusks can be 
grouped into CA and TH ecological groups, as have been 
previously defined in the Chinese Quaternary loess-paleosol 
sequences [49-57]. In the Dongwan section, the TH group 
comprises Gastrocopta armigerella, Gastrocopta sp., Punctum orphana, 
Punctum sp., Metodontia heresowskii, Metodontia huaiensis, Metodontia 
yantaiensis, Metodontia cf huaiensis, Metodontia ci. yantaiensis, Metodontia 
cf heresowskii, Metodontia sp., Kaliella sp., Macrochlamys sp., and Opeas 
sp.; and the CA group comprises Cathaica sp., Cathaica pulveratrix, 
Cathaica pulveraticula, Cathaica schensiensis, Cathaica placenta, Pupilla 
aeoli, Pupilla grahaui, Pupilla sp., Vallonia cf pulchella, Vallonia sp., and 
Pupopsis retrodens. 

The distributions of these terrestrial moUusks in the Dongwan 
section have been described previously [58]. Cathaica sp. is the 
most continuously distributed mollusk taxon in the Dongwan 
section, and Gastrocopta sp., Pupilla sp., and Vallonia sp. are other 
well represented taxa. The TH species of Metodontia, Punctum, 
Macrochlamys, and Opeas are concentrated in the upper part of the 
section. In general, the CA species are dominant in the loess layers, 
while the TH species mainly occur in the paleosols [58]. 

Variations in the number of species (S), diversity [Ft], equitabHity 
(E), and total individuals of all species, CA species, and TH species 
in the Dongwan sequence are plotted against depth in Figure 3. 
These data can be found in the supporting information data 
(Table SI). Variations in these parameters allow definition of four 
moUusk zones. In Zone 1, from the bottom of the sequence to 



about 56.8 m depth (—7.1-6.2 Ma), high values of all the CA 
species (S(CA)) and individuals and diversity (//(CA)) dominate 
over low values of all of the TH species (S(TH)) and individuals 
and diversity (//(TH)). A prominent feature of this zone is that high 
values of S(CA), //(CA), CA mollusk individuals, S(TH) and 
//(TH), as well as total mollusk species and individuals, occur at 
around 70 m depth, at the very base of the sequence; however, 
these values then decrease upwards. //(CA) and //(total) remain at 
high levels throughout the zone, and the equitability (E) of total 
species is also higher than in the other zones. 

In Zone 2, from 56.8-34.8 m depth (-6.2-5.4 Ma), the total 
number of CA species and individuals remains at a similarly high 
level as in the previous zone, except for the interval from 42 to 
38 m. In contrast, the total number of TH species and individuals 
is slightly higher than in the previous zone. The total number of 
species and individuals of all of the species (CA and TH combined) 
are generally higher than in the previous zone but with a 
decreasing trend. The diversity of all species, //(total), does not 
exhibit significant changes compared to Zone 1, and E(total) is 
somewhat lower than in Zone 1. 

Zone 3, from 34.8-16.2 m depth (—5.4-4.4 Ma), is character- 
ized by low numbers of CA species and individuals being almost 
absent in the upper part, from 24 m to 16.2 m. This pattern of 
variation is paralleled by a clear increase of all the TH indices, 
including diversity. The total number of species and individuals of 
all of the species is lower than in Zone 2, except for the interval 
from about 32 to 28 m. There are no large magnitude changes in 
//(total) and E(total) within the zone. 

Zone 4 corresponds to the depth from 16.2 m to the top of the 
sequence (—4.4 to 3.5 Ma). The number of CA species and 
individuals increases markedly and maintains relatively high values 
throughout the zone. The TH indices exhibit high values from 
about 16.2 to 9 m and then decrease significantly to very low 
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values up to 4 m, and increase again above about 4 m depth. All 
of the indices for all species, i.e., S(total), //(total), and E(total), 
increase markedly; however, //(total) and E(total) exhibit a 
generally declining trend, implying that the diversity of the 
terrestrial moUusk populations decreased and that the distribution 
of individuals of different species was uneven. The total moUusk 
individuals in this zone is high and remains relatively stable except 
for the interval between 8 and 4 m, a pattern which differs from 
the other three zones where the total number of moUusk 
individuals is high at the beginning of the zone and then declines 
thereafter. 

Moreover, the variation in diversity of the CA and TH 
ecological groups from the entire Dongwan sequence, as exhibited 
in Figure 3, can be differentiated into two major intervals based 
solely on the Shannon index of the CA and TH groups; and this 
result probably reflects different ecological population dynamics. 
First, from the base of the Dongwan sequence to 34.8 m depth, the 
CA group is dominant, exhibiting high diversity with H values 
varying between 0 and 1.75 (mean of 0.52). In contrast, the 
diversity of the TH populations exhibits significantly low values, 
ranging from 0 to 1.58 (mean of 0.18) within this depth interval, 
although the diversity of all moUusk species, i.e. the sum of CA and 
TH, does not exhibit any clear changes. Second, from 34.8 m 
depth to the top of the sequence, the diversity of the CA group 
declines significantiy with //values varying between 0 and 1.56 
(mean of 0.14); and in contrast the TH group becomes dominant 
in terms of diversity with //values ranging from 0 to 2.64 (mean of 
0.64). The diversity and equitabUity of all moUusk species also 
exhibit a significant change at about 16.2 m depth within this 
interval. 

Discussion 

Late Miocene and Pliocene Paleoclimatic Evolution in the 
Western CLP 

Previous studies of European and North American Quaternary 
terrestrial moUusks have shown that moUusk assemblages, based 
not only on the occurrence of characteristic species but also on the 
statistical dispersion of the assemblages (e.g., diversity, as described 
by the Shannon index (//)), can provide information about climatic 
concUtions and moUusk populations [65-67]. For example, 
correlation of terrestrial moUusk groups with changes in dust flux, 
and thus climatic conditions, has previously been observed in the 
European Qiiaternary terrestrial moUusk diversity record [65]. 
Here we use the Dongwan terrestrial moUusk record to extend the 
statistical analysis of terrestrial moUusk populations and the 
interpretation of climatic changes to the Late Neogene. As shown 
in Figure 4, four stages can be recognized according to changes in 
the moUusk indices used in the present study, and which outUne 
the evolutionary history of ecological and climatic conditions in 
the CLP from 7.1 to 3.5 Ma. 

First, from 7. 1 to 6.2 Ma (Zone 1), the total number of CA 
species (S(CA)) and diversity (//(CA)) index are high; however, the 
number of CA moUusk individuals is generaUy low apart from a 
large high peak at around 7 Ma and which indicates diversified 
CA species. Conversely, aU of the indices of the TH species exhibit 
low values except for the peak at around 7 Ma, indicating low 
diversity. These features may be related to the occurrence of very 
cold, drv' climatic conditions, the occurrence of which is roughly 
coincident with the expansion of C4 vegetation in the northeastern 
CLP [27], Central Inner Mongolia (Tunggur Area, Xilinhot Area, 
and Huade Area) [68] (Figure 1), and Pakistan [2]. They also 
indicate the occurrence of seasons with water stress for vegetation 
and terrestrial moUusks as weU as the fact that relatively dry 



climatic conditions had already appeared in northern China 
during the Late Miocene, coincident in age with the global 
transition from C3 to C4 vegetation [4]. The occurrence of C4 
plants during the latest Miocene should have been restricted to 
limited areas within deserts, and the occurrence of these niches 
could be coincident with, or occurred after, major uplift of the 
Tibetan Plateau at about 8 Ma [3,8,10,11]. However, C4 
vegetation was not yet to develop in the central and southern 
CLP, demonstrating that C3 plants were still dominant as 
indicated by S' ^C \ alues of soU carbonate in the Lingtai, Xifeng 
and Liantian Red Clay sequences in the eastern CLP [6,9,24]. 
This is also evidenced by the thickness of loess layers in the studied 
Dongwan sequence and by dust mass accumulation rates in the 
western Pacific. High dust deposition rates with a graduaUy 
increasing trend culminated at about 6.2 Ma in the western Pacific 
[69], corresponding to thicker loess layers in the Dongwan 
sequence (Figure 4). Both features indicate drier climatic 
conditions in the Asian interior, the potential source region of 
dust deposits in the CLP and western Pacific. WhUe dust 
accumulation in the western Pacific reached a maximum from 
about 6.5 to 6.2 Ma, the number of moUusk species and 
individuals at Dongwan was very low, implying very dry climatic 
conditions which were unsuitable even for the development of CA 
species. Pedostratigraphy and iron geochemistry of the Lingtai 
Red Clay deposits in the CLP indicate that the EA summer 
monsoon was relatively weak from 7.05 to 6.2 Ma [20,21]. The 
5'^C values of fossfl enamel from a diverse group of herbivores 
and of paleosol carbonate and organic matter from the Linxia 
Basin incUcate that C4 grasses were either absent or insignificant in 
the Linxia Basin prior to ~2-3 Ma, suggesting that the East Asian 
summer monsoon was probably not strong enough to affect this 
part of China throughout much of the Neogene [28] . 

However, during our studied interval there is much evidence for 
a strong summer monsoon, evidenced mainly from mammal 
assemblages as well as analyses of mammal tooth and soU 
carbonate isotopes [22-25,27]. Hypsodonty analysis of fossU 
mammals indicates that northern China became more humid at 
7-8 Ma, coincident with the previously recognized onset of Red 
Clay deposition in the eastern CLP, and which was interpreted by 
the authors as reflecting the onset or intensification of summer 
monsoonal precipitation [22,23]. Soil carbonate isotopes from the 
Lantian fluvial and Red Clay deposits indicate the absence of any 
marked change in plant photosynthetic pathway or climate, 
implying the occurrence of pure C3 vegetation with no indications 
of any C4 plants during the Late Miocene and Pliocene. A general 
rise in 8"^0 values probably reflects increased summer precipi- 
tation related to the onset and/or intensification of the Asian 
monsoon system [24,25], which is supported by a mammalian 
faunal turnover event implying a marked change to more humid 
and closed habitats [26]. A strong summer monsoon at 7.1— 
6.2 Ma is also supported by the records of clay/feldspar ratio, 
kaolinite/ [:hlorite ratio and biogenic opal mass accumulation rates 
(MAR) from the SCS, although the authors emphasize that their 
study is a schematic reconstruction which only outiines several 
principal stages and neglects the details [33]. 

There is also much evidence for a weakening of the summer 
monsoon and strengthening of the winter monsoon during the 
Late Miocene. The EA summer monsoon weakened after 7.5 Ma 
as evidenced by combined planktonic foraminifera Mg/Ca and 
5"^0 records from Ocean DriUing Program (ODP) Site 1 146, 
northern SCS [30]. Radiolarian species numbers and individuals 
and diversity from ODP Site 1143, southern SCS, suggest a 
summer monsoon maximum at 8.24 Ma and a major decline after 
7.70 Ma [32]. Geochemical records from die ODP Site 1148, 
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Figure 4. Variations in tlie terrestrial moiluslts of tiie Dongwan sequence in the western CLP from 7.1 to 3.5 Ma, and comparison 
witli other climate proxies. The proxies are of soil carbonate in the Lingtai Red Clay sequence [6], loess and soil thickness in the 

Dongwan sequence (this study), global cooling trend deduced from the marine benthic foraminiferal 5^^0 record [1], and continental aridity inferred 
from the mass accumulation rate (MAR) of eolian dust at ODP Site 885/886 in the western Pacific [69]. Geomagnetic polarity chronology and age of 
the Dongwan sequence is from Hao and Guo (2004) [41]. S(CA)-Total species of the cold-aridiphilous (CA) mollusk group. H(CA)-Diversity of the CA 
mollusk group. S(TH)-Total number of species of the thermo-humidiphilous (TH) mollusk group. H(TH)-Diversity of the TH mollusk group. 
doi:1 0.1 371 /journal.pone.0095754.g004 



northern SCS, indicate tliat chemical weathering intensity 
decreased after the early Miocene with a rapid decrease centering 
at 7.2 IVIa, compatible with a weakening of the EA summer 
monsoon [3 1] . The trend of increasing black carbon concentration 
and accumulation rate at ODP Site 1 1 48 suggests the progressive 
development of the East Asian winter monsoon after about 7 Ma 
[29]. In addition, grain-size evidence from the Linxia Basin of the 
CLP indicates that the EA winter monsoon intensified after 
7.4 Ma and again at 5.3 Ma [19]. 

Thus it is noteworthy that the TH moUusks at Dongwan do 
exhibit peak values at around 7 Ma, suggesting warm, humid 
climatic conditions consistent with a mammalian faunal turnover 
event [26]; however, the peak does not extend to the subsequent 
period, implying that the interval of relatively warm and humid 
climate may have been a relatively brief event within the context of 
an overall very cold and arid climate. In addition, a recent 
synthesis of isotopic and sedimentological analyses, climate 
modeling and an extensive mesowear analysis of the Baode Red 
Clay sequence indicates that the climate at 7.5 Ma was more 
humid than during the youngest interval at 5.7 Ma and that 
variable climatic conditions occurred at ~6.5 Ma. Thus a 
significant decrease in the EA summer monsoon strength from 
7-5.7 Ma in the Baode region, based on three samples dated at 
7.0, 6.5, and 5.7 Ma [70], does not contradict the mollusk results. 
We suggest that much more definitive results would have been 
obtained if a higher sampling density had been used. In addition, 
numerous paleoclimatologists regard the occurrence of Red Clay 
sequences as reflecting aridity or desertification and an intensified 
winter monsoon [16,37,42,43,69,71], which contradicts the 
mammal and isotope record. 

We suggest four possible reasons for the discrepancies amongst 
the extant research results. First, the sensitivity of the proxies used 
may be one of the main reasons; and indeed, it has long been 
found in Qiiaternary loess studies that different climatic proxies 
can respond differently to climatic changes [55]. Second, as 
summarized by Wang et al. (2005) [72] and Kaakinen et al. (2006) 



[24], interpretations of expansions of C4 vegetation are conten- 
tious in terms of whether or not they are a summer monsoon 
proxy [2,4,45]. Third, regional difierence in patterns of climate 
change may have occurred, as suggested by Passy et al. (2009) 
[27]. Finally, low stratigraphic resolution may have resulted in the 
failure to resolve environmental changes in sufficient detail, and 
high resolution mammalian faunal studies (for example using a 20- 
cm interval) may shed light on the observed differences. In 
summary, the conflicts cannot be resolved by the present study and 
future studies employing more sensitive proxies, higher resolution 
or continuous sampling combined with accurate dating may be 
necessary. 

Second, from 6.2 to 5.4 Ma (Zone 2), the dominant CA 
terrestrial moUusks indicate that the climate during this interval 
remained cold and dry, but was not drier than in Zone 1, as 
evidenced by the slightiy greater abundance of TH species and 
individuals. At the beginning of this zone, from about 6.2 to 
5.8 Ma, the number of TH mollusk individuals was higher, and 
this corresponds to a declining rate of dust deposition as indicated 
by the thinness of the loess layers in the Dongwan sequence 
(Figure 4). This supports the finding that decreasing or increasing 
dust deposition affects the development of terrestrial moUusks by 
creating either more- or less- favorable environmental conditions 
[65]. After about 5.8 Ma the number of TH mollusk individuals 
was reduced, which is probably related to the increase of dust 
deposition in the CLP, as shown by the increased thickness of the 
loess represented by this interval (Figure 4). These cold, dry 
climatic conditions were not solely restricted to the western CLP; 
in the eastern CLP, moUusks in the Xifeng Red Clay sequence 
exhibit dominant percentages of the CA group associated with a 
few meso-xerophilous components and a lesser occurrence of TH 
species, also indicating a cold, dry climate [46]. This is also 
reflected by the coarse grain size of the Red Clay sequence in the 
CLP [16], and by the fact that dust deposition rates in the western 
Pacific remained high, albeit lower than before, and with a 
declining trend towards the subsequent time interval from 5.4 and 
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4.4 Ma [69] . The thickness of the Dongwan loess layers was not 
higher than before, except at about 6.2 and 5.6 Ma; however, the 
thickness of the soil layers does not change significandy. The 
values of the soil carbonates in the Lingtai and Xifeng Red Clay 
sequences were similar to those of the previous time interval from 
7.1 to 6.2 Ma, indicating again that C4 plants were not the 
dominant vegetation in the studied area of the CLP [6,9] . These 
Unes of evidence indicate generally cold, dry climatic conditions in 
the EA continent but which were not drier than the lower stage. 
However, the global mean bcnthic foraminifi'ral 8"*0 record 
demonstrates that gloljal ice volume was increasing with the 
highest values occurring at about 5.8-5.7 Ma, followed by a 
decrease towards the upper stage [1]. 

There is much evidence supporting the conclusion from the 
Dongwan moUusk results that the winter monsoon may still have 
been strong from 6.2-5.4 Ma but that the summer monsoon 
strengthens than before although may be to a lesser extent, as 
confirmed by the pedostratigraphy and iron geochemistry of the 
Lingtai Red Clay deposits in the CLP [20,21]. The coeval pollen 
record from the Lingtai section exhibits a predominance of 
Chenopodiaceae and Artemisia, indicating a desert or desert- 
grassland landscape during the Late Miocene [7 3] . A higher dust 
accumulation rate and coarser eolian grain size in the CLP and in 
the North Pacific suggest stronger continental aridity in the Asian 
interior from ~6.2 to ~5 Ma [14,16,69]. The variation in the U- 
ratio of the grain size of the Red Clay, reflecting the changing 
strength of the winter monsoon, indicates a strong winter monsoon 
between 6.1 and 5.4 Ma [13]. The sedimentation rate across the 
CLP, including at the Baode Red Clay sequence, indicates that the 
EA winter monsoon strengthened between 6.26 and 5.4—5.25 Ma 
[14,15]. 

Third, from 5.4 to 4.4 Ma (Zone 3), the prominent feature is 
that all of the CA species and individuals are very few in number 
and in fact they decrease significandy with almost no occurrence 
after about 5 Ma. In contrast, there are high values of the TH 
species numbers and individuals with the maxima occurring at 
around 5-5. 1 Ma when a thick soil layer formed, indicating warm, 
moist climatic conditions. In the eastern CLP, the moUusk fauna 
was characterized by maximum abundance of TH species, an 
absence of CA species, and relatively abundant meso-xerophilous 
taxa [46]. The Xifeng poUen record indicates a significant increase 
of temperate forest plants, also implying humid regional climatic 
conditions during this period [74]. The 5'*C values indicate a 
slight decrease and thus a shift towards more C3 plants in the CLP 
in agreement with the pollen results (Figure 4). Ding et al. (1999) 
identify extremely mature soils in the 5.5-3.85 Ma interval at 
Lingtai, and interpret these as indicating warm and wet climates, 
and possibly the strongest summer monsoons in the past 7 Ma 
[20]. Minimum grain size and sedimentation rate in the CLP 
[15,16], decreased thickness of loess layers in the Dongwan 
sequence (Figure 4), and very low dust deposition rates in the 
western Pacific [69] together suggest \'ery \\arm, moist climatic 
conditions in the CLP and less dry conditions in the Asian interior, 
the potential source area of dust in the CLP and western Pacific. 
However, this warm climate is not clearly recorded by the global 
mean benthic foraminiferal 5'^0 record compiled by Zachos et al. 
(2001) [1]. 

Last, from 4.4 to 3.5 Ma (Zone 4), all of the CA indices increase 
significandy and remain almost unchanged, while in contrast all of 
the TH indices exhibit a different pattern of variabUity with high 
values from 4.4 to 4.0 Ma, coincident with changes in the CA 
indices, an abrupt decrease from 4.0 Ma and increased values 
again from 3.7 Ma. The CA species during this time interval were 
different from those between 7.1 and 5.4 Ma. Indeed, the CA 



species of Cathaica pulveraticula, Cathaica schensiensis, and Pupopsis 
retrodms have only been identified in this zone, as shown by Li et al. 
(2()06a) [58], probably indicating that climatic conditions became 
drier. The coeval Xifeng land snail record from the eastern CLP is 
dominated by meso-xerophUous taxa, associated with a signifi- 
candy reduced abundance of TH species and the paucity of CA 
taxa during the late period [46] . Pollen data from the Xifeng Red 
Clay sequence indicate a typical steppe ecosystem during this 
period [74]. These variations are in good agreement with the 
expansion of C4 plants in the central CLP as observed at the 
Lingtai and Xifeng Red Clay sequences [6,9]. In addition, the 
loess layers in the Dongwan section during this time interval are of 
moderate thickness, thinner than from 7.1 to 6.2 Ma but thicker 
than from 5.4 to 4.4 Ma. In contrast the paleosols are somewhat 
thicker than before (Figure 4), indicating longer pedogenesis under 
different climate conditions, and probably enhanced seasonality, 
corresponding to expansions of C4 plants in the central CLP [6,9]. 
Pedostratigraphy and iron geochemistry of the Lingtai Red Clay 
deposits in the CLP indicate that the EA summer monsoon 
weakened significandy over the interval 3.85 to 3.15 Ma [20,21]. 
Dust deposition increased again in the Pacific in parallel with 
increased aridity in the dust source areas [69] (Figure 4). 

Moreover, the entire Dongwan terrestrial moUusk record 
provides information about the evolution of moUusk diversity 
and climatic changes from the Late Miocene to the Pliocene. 
Previous studies indicate that in most Quaternary loess sequences 
in Europe and North America generally fewer than 20 moUusk 
species can be identified and diversity varies with H values 
between 0 and 4, suggesting that the loess deposits provided few 
additional ecological niches for land snaUs to grow and develop 
[65-67]. In the Neogene Dongwan loess sequence in the CLP, 
East Asia, 24 moUusk species were identified and diversity varies 
between 0 and 2.5, which also suggests that fewer niches were 
provided. The response of the Neogene terrestrial moUusk 
populations to climatic changes depends on their ecological 
requirements. Different ecological groups, such as CA and TH, 
respond differently. During the time interval from 7.1 to 5.4 Ma 
when a cold, dry climate obtained, as indicated by a high flux of 
dust in the CLP and western Pacific, the diversity of the CA group 
was high, indicating that relatively cold, dry climatic conditions 
may be favorable for the development of CA terrestrial moUusk 
populations. In contrast, during the r(;lati\'ely warm, moist time 
interval from 5.4 to 4,4 Ma, corresponding to a reduced loess 
thickness in the CLP and low dust flux in the western Pacific, the 
diversity of the TH terrestrial moUusk populations was high. 

Possible Causes of Paleoecological and Paleoclimatic 
Evolution in the Western CLP during the Late Miocene 
and Pliocene 

Dust in the Asian interior, including northwestern China, was 
emitted and transported in two possible modes, corresponding to 
patterns of low and high level atmospheric circulation. One mode 
is that the Asian dust was transported eastwards by high 
atmospheric circulations (westerlies) and reached northwestern 
Pacific, as recorded by high values of dust flux to the northwestern 
Pacific. The other mode is that dust in the Asia interior was 
transported by low \v\cX atmospheric circulation, the EA winter 
monsoon, to the middle reaches of the Yellow River, leading to the 
formation of the CLP [37,49,75,76]. Thus loess sequences in the 
CLP and dust deposits in the western Pacific both relate to climatic 
changes that impacted the regions to the north of the Tibetan 
Plateau, i.e., the Asian interior, as has been indicated by the 
previous studies of Hovan et al., (1989) [77] and Rea et al, (1998) 
[69]. 
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At approximately 8-7 Ma, the accumulation rate of eolian 
deposits in the CLP reached high values of about 4 cm/ka [37] 
and dust deposition in the western Pacific was maximal [69] 
(Figure 4), indicating that particularly dry climatic conditions must 
have prevailed in the Asian interior resulting in the mobilization 
and transport of large amounts of dust. Indeed, a palynological 
study of the Late Miocene-Pliocene sediments of the Dushanzi 
section from northwestern China (Figure 1) indicates that steppe 
taxa [Artemisia and Chenopodiaceae) were generally dominant, 
implying that a dry climate existed in the inland basins of 
northwestern China since 8.7 Ma, except for a warm and humid 
phase that lasted from 5.8 to 3.9 Ma [78]. In addition, significant 
desertification in northwestern China prevailed as early as 7.2- 
7 Ma, as shown by the development of eolian sand dunes in the 
Taklimakan Desert [79] and eolian Red Clay deposition around 
the Lanzhou region [80]. These climatic and environmental 
conditions probably provided a suitable environment for the 
growth and development of C4 plants in northern China and 
promoted the expansion of eolian deposits from the west towards 
the eastern part of the CLP [42,43] and which constituted suitable 
environments for terrestrial moUusks. Thus terrestrial moUusks 
grew and developed in the CLP and were able to record the 
patterns of environmental evolution as recorded by the Dongwan 
sequence. 

Variations in the Dongwan terrestrial moUusk assemblages 
indicate that, during the time interval from 7.1 to 3.5 Ma, major 
paleoecological and paleoclimatic changes occurred at about 6.2, 
5.4, and 4.4 Ma. The change at around 6.2 Ma, the boundary 
betvv("i'n Zones 1 and 2, is less significant than the one at 4.4 Ma, 
and therefore it is not clearly reflected in terms of a transition in 
periodicity recorded by the relative abundance of moUusks at this 
time [47]. Indeed, a relatively high number of CA moUusk species 
and individuals, consistent with the high relative abundance of C A 
moUusks [47], prevailed during Zone 1 and Zone 2, indicating that 
a cold, arid climate prevaUed in the study area. These two zones 
exhibit the same dominant 100 kyr periodicity [47], again 
suggesting that the 6.2 Ma datum was not particularly significant 
at Dongwan. However, the TH moUusk individuals and species 
are somewhat different within these two zones, and in addition the 
a])pc'aran( e of moUusks in the Xifeng Red Clay sequence and the 
upper age of the QA-I section imply that 6.2 Ma may be an 
important datum in the paleoclimatic evolution of the CLP 
[37,46]. 

Changes in the CA moUusks of Zones 1 and 2 are paralleled by 
a global cooling trend [1] (Figure 4), increased ice-rafted detrital 
flux in the Northern Hemisphere [81-85], and the buUdup of the 
Western Antarctic ice sheet [1]. The CLP is particularly sensitive 
to changes in high northern latitudes through the EA winter 
monsoon circulation [76,86] . Extended ice sheets in the Northern 
Hemisphere reinforce the southward movement of cold air and 
thereby enhance the Siberian High that controls the EA winter 
monsoon wind system [86]. Thus, global cooling, especially the 
expansion of ice deposits in high northern latitudes, could have 
affected the Siberian High and the RA winter monsoon, thereby 
expanding habitats for the CA moUusks in the CLP. 

There is no evidence for major uplift of the Tibetan Plateau at 
about 6.2 Ma, and therefore this mechanism is excluded as a cause 
of the change at this time. However, if the Tibetan Plateau 
reached a significant height at about 8 Ma it would have 
accelerated climatic cooling and strengthened the EA winter 
monsoon [3,8,87-90], causing cold, dry climatic conditions in the 
western CLP during the Late Miocene, as discussed above. 

The shift at about 5.3-5.4 Ma roughly corresponds to the onset 
of the PUocene when global climate became warmer than before. 



However, the climatic drivers contributing to the Pliocene global 
warming are stiU highly debated. As summarized by Haywood et 
al. (2009) [91], possible candidates include paleogeographic 
changes [92], altered atmospheric trace gas concentrations and 
water vapor content [93], changes in oceanic circulation [94,95], 
oceanic heat transport [5,96], thermal structure of the oceans [97— 
99], and feedbacks generated through altered land cover 
(including ice sheet extent), surface albedo, cloud cover and 
temperature [100]. 

Our results seem to support the suggestion that the warming 
interval may be related to changes in ocean circulation and ocean 
heat transport caused by the closures of the Panama and 
Indonesian seaways [5,94—96]. Closure of these two seaways 
may have caused changes in heat distribution between the Pacific 
and Atiantic, causing reorganization of global climatic patterns 
and changing the pattern of atmospheric moisture flux from 
latitudinal to meridional, resulting in increased moisture flux to 
high latitudes [5,95,101] and contributing to climate changes in 
the CLP. Moreover, both geological records and modeling studies 
show that closures of the Panama and Indonesian seaways likely 
played important roles in the strengthening and enlargement of 
the western Pacific warm pool [5,102-106], providing more 
moisture and heat to the CLP favorable for the abundant 
occurrence of TH moUusk species. 

The 4.4 Ma shift, representing a climatic transition from early 
Pliocene warming to late Pliocene cooling, is obser\'ed in the 
Dongwan moUusk diversity and species numbers and individuals; 
however, it is not clearly shown in the relative abundance of 
moUusks [47] . It may be related to tectonic events such as the uplift 
of Tibetan Plateau, since climate models suggest that uplift played 
a particularly important role in the evolution of the global 
paleoenvironment [8,87-90], and the effects of which on the EA 
winter monsoon are more significant than on the summer 
monsoon [89]. Changes in depositional facies from distal aUuvial 
plains to proximal alluvial fans and an increase in sedimentation 
rate near Yecheng (Figure 1) in the western Kunlun Mountains 
indicate the uplift of the northern Tibetan Plateau at about 
4.5 Ma [7] . This tectonic activity is thought to trigger not only the 
enhancement of the EA monsoon, but is also considered to be a 
driver of general cooling through the consequent increase in the 
rate of chemical weathering, thus accelerating ice expansion in 
Northern Hemisphere high latitudes [88,107,108] and further 
strengthening the Siberian High and the EA winter monsoon and 
transportation of dust to the CLP, Furthermore, uplift of the 
Tibetan Plateau could have blocked moisture transport to the 
Asian interior, contributing to its aridification. Under these 
climatic conditions, the northern CLP would have commenced a 
drying trend earUer than the southern CLP, causing significant 
ecological changes in terrestrial ecosystems including terrestrial 
moUusks and C4 plants at roughly 4.4 Ma. 

Conclusions 

The Dongwan terrestrial moUusk record from the western CLP 
exhibits four stages during the time interval from 7.1 to 3.5 Ma, 
indicating the phased evolution of paleoecology and paleoclimate. 
From 7.1 to 6.2 Ma, very cold and dry climatic conditions 
prevaUed. From 6.2 to 5.4 Ma, the climate remained cold and dry 
but was not as dry as the preceding interval, as evidenced by the 
dominance of CA moUusks and rather more TH species and 
individuals. From 5.4 to 4.4 Ma, very warm and moist climatic 
conditions prevailed, evidenced by high values of the TH species 
and individuals, as weU as by the very smaU numbers of CA species 
and individuals and their almost complete absence after about 
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5 Ma. From 4.4 to 3.5 Ma, all of the CA indices increase 
significantly and remain at a high level, indicating a cooling 
climate; however, aU of the TH indices exhibit relatively high 
valuesfrom 4.4 to 4.0 Ma, an abrupt decrease from 4.0 Ma and an 
increase again from 3.7 Ma. The three CA species, Cathaka 
pulveraticula, Cathaica schensiensis, Pupopsis retrodens, occurred solely 
during this period and were absent from 7.1 to 5.4 Ma, suggesting 
that the climate from 4.4 to 3.5 Ma was becoming colder and drier 
than previous stage. 

The very cold and arid climatic conditions, with changes at 
about 6.2 Ma, are paralleled by a global cooling trc'nd, increased 
in ice-rafted detrital flux in the Northern Hemisphere [81-85], 
and the buildup of the Western Antarctica ice sheet [1]. The shift 
at about 5.3-5.4 Ma roughly corresponds to the onset of the 
Pliocene when global climate became increasingly warmer than 
previously. Our results si'(^m to support the suggestion that the 
warming interval may be related to changes in oceanic circulation 
and oceanic heat transport [5,94—96]. The 4.4 Ma shift may be 
related to tectonic events such as the uplift of Tibetan Plateau. 

Variations in the diversity of the GA and TH moUusks, //(CA) 
and //(TH), are closely related to chmatic changes during the Late 
Miocene to Pliocene. From 7.1 to 5.4 Ma when a cold, dry climate 
prevailed, the CA group was the more diverse. In contrast, during 
the relatively warm, wet time interval from 5.4 to 4.4 Ma, the TH 
terrestrial moUusk populations became more diverse. It should be 
pointed out that most of the Neogene terrestrial moUusks in the 
CLP have modern representatives and therefore they have the 
potential to estimate quantitatively Neogene changes in temper- 
ature and precipitation in the CLP. However, such estimates 
depend upon the development of a training set Ijascd on a large 
number of surface samples. Changes in fossil terrestrial moUusk 
diversity are significant for the prediction of terrestrial biodiversity 
changes, and Quaternary loess deposits in Europe and North 
America have been studied in this context. The 22-Ma loess 
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deposits in China provide an excellent opportunity for under- 
standing long term changes in terrestrial mollusk diversity; and 
although the present paper focuses on the Late Neogene terrestrial 
mollusk in the CLP, ongoing studies will focus on other time 
intervals. 

Supporting Information 

Table SI The depth, age, and mollusk data &om the 
Dongwan loess-paleosol sequence. S(CA) Total number of 

species of the cold-aridiphilous (CA) mollusk group. H{GA)- 
Diversity of the CA mollusk group. CAMI-MoUusk individuals of 
the CA mollusk group. S(TH)-Total number of species of the 
thermo-humidiphilous (TH) moUusk group. //(TH)-Diversity of 
the TH mollusk group. THMI— Mollusk individuals of the TH 
mollusk grotip. S(total)-Total number of species of the total 
mollusk group. //(total)-Diversity of the total moUusk group. 
E(total)-Equitability of the total moUusk group. TMI-Total 
MoUusk individuals. 
PCLS) 

Acknowledgments 

We thank all of the reviewers for valuable comments, criticisms and 

suggestions. We are grateful to Dr. Jan Bloemendal for editing the English 
and for suggestions, Prof. Qingzhen Hao for field assistance and Prof 
Shiling Yang for pro\'iding us with 6^^C data of soil carbonates from the 
Lingtai Red Clay sequence. 

Author Contributions 

C^oiieeived and designed the experiments: FL NW. Performed the 
experiments: FL NW YP. Analyzed the data: FL NW DDR YD DZ. 
Contributed reagents/ materials/ analysis tools: FL NW. Wrote the paper: 
FL NW DDR. 



14. Wen LJ, Lu HY, Qiang XK (2005) Changes in grain-size and sedimentation 
rate of the Neogene Red Clay deposits along the Chinese Loess Plateau and 
implications for the palaeowind system. Sci China Earth Sci 48: 1452-1462. 

15. Zhu YM, Zhou LP, Mo DW, Kaakinen A, Zhang ZQ, et al. (2008) A new 
magnetostratigraphic framework for late Neogene Hipparioii Red C^lay in the 
eastern Loess Plateau of China. Palaeogeogr. Palaeoclimatol. Palaeoecol. 268: 
47-57. 

16. Guo ZT, Peng SZ, Hao QZ, Biscaye PE, An ZS, et al. (2004) Late Miocene- 
Pliocene development of Asian aridification as recorded in the Red-Earth 
Formation in the northern China. Glob. Planet. Change 41: 135—145. 

17. Wang YX, YangJD, Chen J, Zhang KJ, Rao WB (2007) The Sr and Nd 
isotopie variations of the Chinese Loess Plateau during the past 7 Ma: 
Implications for the East Asian winter monsoon and source areas of loess. 
Palaeogeogr. Palaeoclimatol. Palaeoecol. 249: 351-361. 

18. Clhen J, Chen Y, Liu L\V, Ji JF, Balsam W, et al. (2006) Zr/Rb ratio in die 
C^hinese loess secjuences and its implication for changes in the East Asian winter 
monsoon strength. Geochim. Cosnioehini. Acta 70: 1471-1482. 

19. Fan MJ, Song GH, Dettman DL, Fang XM, Xu XH (2006) Intensification of 
the Asian winter monsoon after 7.4 Ma: Grain-size evidence from the Linxia 
Basin, northeastern Tibetan Plateau, 13.1 Ma to 4.3 Ma. Earth Planet. Sci. 
Lett. 248: 186-197. 

20. L:iing ZL. Xiong SF, Sun JM, Yang SL, Gu ZY, et al. (1999) Pedostratigraphy 
and pak:omagnetism of a 7.0 Ma eoliaii loess-red clay sequence at Lingtai, 
Loess Plateau, north-central China and the implications for palcomonsoon 
evolution. Palaeogeogr. Palaeoclimatol. Palaeoecol. 152: 49-66. 

21. Ding ZL, Yang .SL, Sun JM, Liu TS (2001) Iron geochemistry of loess and red 
elav deposits in the Chinese Loess Plateau and implications for long-term Asian 
monsoon evolution in the last 7.0 Ma. Earth Planet. Sci. Lett. 185: 99-109. 

22. Fortelius M, Eronen JT, Jernvall J, Liu L, Pushkina D, et al. (2002) Fossil 
mammals resolve regional patterns of Eurasian climate change during 20 
million years. Evol. Ecol. Res. 4: 1005-1016. 

23. Liu LP, Eronen JT, Fortelius M (2009) Significant Mid-Latitiide Aridity in die 
Middle Miocene of East Asia. Palaeogeogr. Palaeoclimatol. Palaeoecol. 279: 
201-206. 

24. Kaakinen A, Sonninen E, Lunkka JP (2006) Stable isotope record in paleosol 
carbonates from the Chinese Loess Plateau: implications for late Neogene 



PLOS ONE I www.plosone.org 



10 



April 2014 | Volume 9 | Issue 4 | e95754 



Late Miocene-Pliocene Paleoclimate in Chinese Loess Plateau 



palcoclimatc and paleovegetation. Palaeogeogr. PalaeoclimatoL PalaeoecoL 

237: 359-369. 

25. Kaakinen A, Passey BH, Zhang Z, LiuL, Pesonen LJ, et al. (2013) Stratigraphy 
and Paleoecology of the classical dragon bone localities of Baode County, 
Shanxi Province, p. 203-217. In Wang X, et al. (eds): Fossil MammEils of Asia: 
Neogene Biostratigraphy and Chronology. Columbia University Press, New 
York. 

26. Zhang ZQ, (}cntr\ AW, Kaakincn A, Liu LP, Lunkka JP, ct al. (2002) Land 
mammal faunal sequence of" the late Miocene of China: new evidence from 
Lantian, Shaanxi Province. Vert. PalAsiat. 40: 165-176. 

27. Passey BH, Ayliife LK, Ka^nen A, Zhang ZQ, Eronen JT, et al. (2009) 
Strengthened East Asian summer monsoons during a period of high-latitude 
warmth? Isotopic evidence from Mio-Phocene fossil mammals and soil 
carbonates from northern China. Earth Planet. Sci. Lett. 277: 443-^52. 

28. Wang Y, Deng T (2005) A 25 m.y. isotopic record of paleodiet and 
environmental change from fossil mammals and paleosols from the NE margin 
of the Tibetan Plateau. Earth Planet. Sci. Lett. 236: 322-338. 

29. Jia GD. Peng PA, Zhao QH, Jian ZM (2003) Changes in terrestrial ecosystem 
since 30 JVla in East Asia: Stable isotope evidence from black carbon in the 
South China Sea. Geology 31: 1093-1096. 

30. Steinke S, Groeneveld J, Johnstone H, Rendle-Biihring R (2010) East Asian 
summer monsoon weaJcening after 7.5 Ma: Evidence from combined 
planktonic foraminifera Mg/Ca and 6^^0 (ODP Site 1146; northern South 
China Sea). Palaeogeogr. PalaeoclimatoL PalaeoecoL 289: 33^3. 

31. Wei GJ, Li XH, Liu Y, Shao L, Liang X (2006) Gcochcmieal record of 
chemical weathering and monsoon climate change since the early Miocene in 
the South China Sea. Palcoceanography 21: PA4214, doi:10.1029/ 
2006PA001300. 

32. Chen MH, Wang RJ, Yang LH, Han JX, Lu J (2003) Development of east 
Asian summer monsoon environments in the late Miocene: radiolarian 
evidence from Site 1143 of ODP Leg 184. Mar. Geol. 201: 169-177. 

33. Wan SM, Li AC, Clift PD, Jiang H (2006) Development of the East Asian 
summer monsoon: Evidence from the sediment record in the South China Sea 
since 8.5 Ma. Palaeogeogr. PalaeoclimatoL PalaeoecoL 241: 139-159. 

34. Jacques FMB, Shi G, Wang WM (2013) Neogene zonal vegetation of China 
and the evolution of the winter monsoon. Bull. Geosci. 88: 175-193. 

35. Liu YS, L'teseher '1', Zhou ZK, Sun BN (2011) The evolution of Miocene 
climates in North China: Preliminary results of quantitative reconstructions 
from plant fossil records. Palaeogeogr. Palaeocliniatol. PalaeoecoL 304: 308- 
317. 

36. Tang H, Micheels A, Eronen JT, ForteHus M (201 1) A regional climate model 
experiment to investigate the Asian monsoon in the Late Miocene. CHm. Past 

7: 847-868. 

37. Guo ZT, Ruddiman WF, Hao QZ, Wu HB, Qiao YS, et al. (2002) Onset of 
Asian desertification by 22 Myr ago inferred from loess deposits in China. 
Nature 416: 159-163. 

38. Sun XJ, Wang PX (2005) How old is the Asian monsoon system?- 
Palaeobotanical records from China. Palaeogeogr. PalaeoclimatoL PalaeoecoL 

222: 181-222. 

39. Qiang XK, An ZS, Song YG, Chang H, Sun YB, et al. (201 1) New eoHan red 
clay sequence on the western Chinese Loess Plateau linked to onset of Asian 
desertification about 25 Ma ago. Sci China Earth Sci 54: 136-144. 

40. Hao QZ, Guo ZT (2007) Magnetostratigraphy of an early-middle Miocene 
loess-soil sequence in the western Loess Plateau of China, (ieophys. Res. Lett. 
34: Lia305, doi:10.1029/2003C;L031 162. 

41. Hao Q_Z, (juo ZT (2004) Magnetostratigraphy of a late Miocene-Pliocene 
loess-paleosol sequence in the western Loess Plateau in China. Geophys. Res. 
Lett. 31: L09209, doi:10.1029/2003GL019392. 

42. Sun DH, Shaw J, An ZS, Chen MY, Yue LP (1998) Magnetostrat%raphy and 
paleoclimatic interpretation of a continuous 7.2 Ma late Cenozoic eoHan 
sediments from the Chinese Loess Plateau. Geophys. Res. Lett. 25: 85-88. 

43. Ding ZL, SunJM, Yang SL, Liu TS (1998) Preliminary magnetostratigraphy of 
a thick eolian red clay-Loess sequence at Lingtai, the Chinese Loess Plateau. 
Geophys. Res. Lett. 25: 1225-1228. 

44. Xu Y, Yue LP, Li JX, Sun L, Sun B, et al. (2009) An 11 -Ma-old red clay 
sequence on the Eastern Chinese Loess Plateau. Palaeogeogr. PalaeoclimatoL 
PalaeoecoL 284: 383-391. 

45. An ZS, Huang YS, Liu WG, Guo ZT, Clemens SC, et al. (2005) Multiple 
expansions of C4 plant biomass in East Asia since 7 Ma coupled with 
strengthened monsoon circulation. Geology 33: 705-708. 

46. Wu NQ, Pei YP, Lu HY, Guo ZT, Li FJ, et al. (2006) Marked ecological shifts 
during 6.2—2.4 Ma revealed by a terrestrial moUuscan record from the Chinese 
Red Clay Formation and implication lor palacoclimatic evolution, Palaeo- 
geogr. PalaeoclimatoL PalaeoecoL 233: 287-299. 

47. Li IJ, Rousseau DD, Wu NQ, Hao QZ, Pei YP (2008) Late Neogene evolution 
of the East Asian monsoon revealed by terrestrial mollusk record in western 
Chinese Loess Plateau: from winter to summer dominated sub-regime. Earth 
Planet. ScL Lett. 274: 439-^47. 

48. Zoller L, Semmel A (2001) 175 years of loess research in Germany-long 
records and "unconformities". Earth-ScL Rev. 54: 19-28. 

49. Liu TS (1985) Loess and the Environment. China Ocean Press, Beijing, 251 p. 

50. Rousseau DD, Wu NQ_ (1997) A new moUuscan record of the monsoon 
variability over the past 130 000 yr in the Luochuan loess sequence, China. 
Geology 25: 275-278. 



51. Rousseau DD, Wu NQ_(1999) Mollusk record of monsoon variability during 
the L2-S2 cycle in the Luochuan loess sequence, China. Quat. Res. 52: 286- 
292. 

52. Rousseau DD, Wu NQ, Guo ZT (2000) The terrestrial mollusks as new indices 
of the Asian pjileomonsoons in the Chinese loess plateau. Glob. Planet. Change 
26: 199-206. 

53. Wu NQj Rousseau DD, Liu TS (1996) Land mollusk records from the 
Luochuan loess sequence and their paleoenvironmental significance. ScL 

China, Ser. D: Earth Sci. 39: 494-502. 

54. Wu NQ, Rousseau DD, Liu XP (2000) Response of mollusk assemblages from 
the Luochuan loess section to orbital forcing since the last 250 ka. Chin. Sci. 
Bull. 45: 1617-1622. 

55. Wu NQ, Rousseau DD, Liii TS, Lu HY, Gu ZY, et al. (2001) Orbital forcing of 
terrestrial mollusks and elimalie changes from the Loess Plateau of China 
during the past 350 ka. J. Geophys. Res. 106: 20045-20054. 

56. Wu NQ, Liu TS, Liu XP, Gu ZY (2002) Mollusk record of millennial dimate 
variability in the Loess Plateau during the Last Glacial Maximum. Boreas 3 1 : 
20-27. 

57. Wu NQ Chen XY, Rousseau DD, Li IJ, Pei YP, ct al. (2007) Climatic 
conditions recorded by terrestrial moUusc assemblages in the Chinese Loess 
Plateau during marine Oxygen Isotope Stages 12—10. Quat. Sci. Rev. 26: 

1884-1896. 

58. Li FJ, Wu NQ, Pei YP, Hao QZ, Rousseau DD (2006a) Wind-blown origin of 
Dongwan late Miocene-Pliocene dust sequence documented by land snail 
record in western Chinese Loess Plateau. Geology 34: 405^08. 

59. Li FJ, Wu NQ, Rousseau DD (2006b) Preliminary study of mollusk fossils in the 
Qinan Miocene loess-soU sequence in western Chinese Loess Plateau. Sci. 
China, Ser. D: Earth Sci. 49: 724-730. 

60. Hou XY (1983) Vegetation of China with reference to its geographical 
distribution. Annals of the Missouri Botanical Garden 70: 509-549. 

61. Cande SC, Kent DV (1995) Revised calibration of the geomagnetic polarity 
timescale for the Late Cretaceous and Cenozoic. J. Geophys. Res. 100: 6093— 
6095. 

62. Kukla G, An ZS, Melice JL, Gavin J, Xiao JL (1990) Magnetic susceptibility 
record of Chinese loess. Trans. R. Soc. Edinb. Earth Sci. 81: 263-288. 

63. Puissegur JJ (1976) MoUusques continentaux quatemaires de Bourgogne. 
Significations stratigraphiques et climatiques. Rapports avec d'autres faunes 
boreales de France: Universite de Dijon Memoires Geologiques, 3: 241 p. 

64. Shannon (^E (1948) A mathematical theory of communication. Bell Syst. tech. 
J. 27: 379 423. 

65. Rousseau DD (1992) Terrestrial mollusks as indicators of global aeolian dust 
fluxes during glacial stages. Boreas 21: 105—109. 

66. Rousseau DD, Limondin N, Puissegur JJ (1993) Holocene environmental 
signals from mollusk assemblages in Burgundy (France). Quat. Res. 40: 237— 
253. 

67. Rousseau DD, Kukla G (1994) Late Pleistocene climate record in the Eustis 
loess section, Nebraska based on land snail assemblages and magnetic 
susceptibily. Quat. Res. 42: 176-187. 

68. Zhang CF, Wang Y, Deng T, Wang XM, Biasatti D, et al. (2009) C4 expansion 
in the central Inner Mongolia during the latest Miocene and early Pliocene. 
Earth Planet. Sci. Lett. 287: 31 1-319. 

69- Rea DK, Snoeekx H, Joseph LH (1998) Late CVnozoie eolian deposition in the 
North Pacific: Asian drying, Tibetan uplift, and cooling of the northern 
hemisphere. Paleoeeanography 13: 215-224. 

70. Eronen JT, Kaakincn A, Liu LP, Passey BH, Tang H, et al. (2014) Here be 
Dragons : Mesowear and tooth enamel isotopes of the classic Chinese 
"Hipparion" fauuEis from Baode, Shanxi Province, China. Ann. ZooL Fennici 

51: 227-244. 

71. Zheng HB, PoweU CM, Rea DK, WangJL, Wang PX (2004) Late Miocene 
and mid-Pliocene enhancement of the East Asian monsoon as viewed from the 
land and sea. Glob. Planet. Change 41: 147-155. 

72. Wang PX, Clemens S, Beaufbrt L, Braconnotd P, Ganssen G, et al. (2005) 
Evolution and variability of the Asian monsoon system: state of the art and 
outstanding issues. Quat. Sci. Rev. 24: 595-629. 

73. Wu YS (2001) Palynoflora at late Miocene-early Pliocene from Leijiahe of 
Lingtai, Gansu Province, China. Acta Botanica Sinica 43: 750-756 (in Chinese 
with English abstract). 

74. Wang L, Lu HY, Wu NQ, Li J, Pei YP, et al. (2006) Palynological evidence for 
Late Miocene-Pliocene vegetation evolution recorded in the red clay sequence 
of the central Chinese Loess Plateau and implication for palaeoenvironmental 
change. Palaeogeogr. PalaeoclimatoL PalaeoecoL 211: 118 128. 

75. An ZS, Liu TS, Lu Y(J, Porter SC, Kukla (J. el al. il990l The long-term 
paleomonsoon variation recorded by the loess-paleosol sequence in Central 
China. Quat. Int. 7: 91-96. 

76. Ding ZL, Liu TS, Rutter NW, Yu ZW, Guo ZT, et al. (1995) Ice-volume 
forcing of East Asian winter monsoon variations in the past 800,000 years. 
Quat. Res. 44: 149-159. 

77. Hovan SA, Rea DK, Pisias NG, Shackleton NJ (1989) A direct hnk between die 
China loess and marine 5'^0 records: Aeolian flux to the North Pacific. Nature 
340: 296-298. 

78. SunJM, Xu QH, Huang BC (2007) Late Cenozoic magnetochronology and 
paleoenvironmental changes in the northern foreland basin of the Tian Shan 
Mountains. J. Geophys. Res. 112: B04107, doi:10.1029/2006JB004653. 



PLOS ONE I www.plosone.org 



11 



April 2014 | Volume 9 | Issue 4 | e95754 



Late Miocene-Pliocene Paleocllmate in Chinese Loess Plateau 



79. Sun JM, Zhang ZQ, Zhang LY (2009) New evidence on the age of the 
Taklimakan Desert. Geology 37; 159—162. 

80. Sun DH, Zhang YB, Han F, Zhang Y, Yi ZY, et al. (2011) Magnetostrati- 
graphy and palaeoenvironmental records for a Late Cenozoic sedimentary 
sequence from Lanzhou, Northeastern margin of the Tibetan Plateau. Glob. 
Planet. Change 76: 106-116. 

81. Janscn E, SjoholmJ (1991) Reconstruction of glaciation over the past 6 Myr 
from icc-bornc deposits in the Norwegian Sea. Nature 349: 600—603. 

82. dcMenoeal PB (1993) Wireline logging on the North Pacific transect. JOIDES 
J. 19: 29. 

83. Wolf-Welling TCW, Gremer M, O'ConncU S, WinUcr A, Thiede J (1996) 
Cenozoic Arctic gateway paleocllmate variability: indications from changes in 
coarse-fraction composition. Proc. Ocean Drill. Program, Sci. Results 151: 

515-567. 

84. Thiede J, Winlder A, Wolf-Weilling T, Eldholm O, Myhre AM, et al. (1998) 
Late Cenozoic history of the polar north Atlantic: results irom ocean drilling. 
Qiiat. Sci. Rev. 17: 185-208. 

85. St. John KEK, Krissek L.A (2002) The late Miocene to Pleistocene ice-rafting 
historv of southeast (Greenland. Boreas 31: 28-35. 

86. Hao QZ, Wang L, Oldiicld F, Peng SZ, Qin L, et al. (2012) Delayed build-up 
of Arctic ice sheets during 400000-year minima in insolation variability. Nature 
490: 393-396. 

87. Kutzbaeh JE, Guetter PJ, Ruddiman WF, Prell WL (1989) Sensitivity of 
climau' lo late cenozoic uplift in southern Asia and American West: Numerical 
experiments. J. Geophys. Res. 94: 18393-18407. 

88. Ruddiman WF, KutzbachJE (1990) Late Cenozoic plateau uplift and climate 
change. Trans. R. Soe. Edinb. Earth Sei. 81: 301-314. 

89. Liu XD, Yin ZY (2002) Sensitivity of East Asian monsoon climate to the uplift 
of the Tibetan Plateau. Palaeogeogr. Palaeoclimatol. Palaeoecol. 183: 223-245. 

90. Liu XD, KutzbachJE, Liu ZY, An ZS, Li L (2003) The Tibetan Plateau as 
amplifier of orbital-scale variability of the East Asian monsoon. Geophys. Res. 
Lett. 30: 16, 1839, doi:10.1029/2003GL017510. 

91. Ha>™ood AM, Dowsett HJ, Valdes PJ, Lunt DJ, Francis JE, et al. (2009) 
Introduction. Pliocene climate, processes and problems. Phil. 'Frans. R. Soe. A, 
367: 3-17. 

92. Rind D, Chandler M (1991) Increased ocean heat transports and warmer 
climate. J. Geophys. Res. 96: 7437-7461. 

93. Raymo ME, Rau G.H (1992) Plio-Pleistocene atmospheric CO2 levels inferred 
from POM 8'^C at DSDP Site 607. Eos Trans. AGU 73 (SuppL), 95. 

94. Ravelo AC, Andreasen DH (2000) Enhanced circulation during a warm 
period. Geophys. Res. Lett. 27: 1001-1004. 



95. Cane MA, Molnar P (2001) Closing of the Indonesian seawav as a precursor to 
cast African aridifieation around 3—4 million years ago. Nature 411: 157—162. 

96. Dowsett HJ, Cronin TM, Poore PZ, Thompson RS, Whadey RC, et al. (1992) 
Micropaleontological evidence for increased meridional heat-transport in the 
North Adantic Ocean dtu-ing the Pliocene. Science 258: 1133-1135. 

97. Philander GS, Fedorov AV (2003) Role of tiopics in changing the response to 
Milankovich forcing some three million years ago. Paleoceanography 18: 
PA1045, doi:10.1029/2002PA000837. 

98. Wara MW, Ravelo AC, Delaney ML (2005) Permanent El Nino-like conditions 
during the Pliocene warm period. Science 309: 758—761. 

99. Fedorov AV, Dckens PS, McCarthy M, Ravelo AC, deMenocal PB, et al. 
(2006) The Pliocene paradox (mechanisms for a permanent El Nino). Science 
312: 148.5-1489. 

100. Haywood A.\l, Valdes PJ (2004) Modelling Middle PUoeene warmth: 
contribution of atmosphere, oceans and eryosphere. Earth Planet. Sci. Lett. 
218: 363-377. 

101. Hall R (2002) Cenozoic geological and plate tectonic evolution of SE Asia and 
the SW Pacific: Computer-based reconstructions and animations. J. Asian 
Earth Sci. 20: 353-434. 

102. Maier-Reimer E, Mikolajewiez U, Crowley T (1990) Ocean general circulation 
model sensitivity experiment with an open central American isthmus. 
Paleoceanography 5: 349—366. 

103. Mikolajewiez U, Crowley IJ (1997) Response of a coupled ocean/energy 
balance model to restricted How through the central American isthmus. 
Paleoceanography 12: 429-441. 

104. Mikolajewiez U, Maier-Reimer E, Crowley TJ, Kim KY (1993) Effect of Drake 
and Panamanian gateways on the circulation of an ocean model. Paleoceano- 
graphy 8: 409-^26. 

105. Chaisson WP, Ravelo AC (2000) Pliocene development of east-west 
hydrographic gradient in the equatorial Pacific. Paleoceanography 15: 497- 
505. 

106. Li BH, Wang JL, Huang BQ, Li QY, Jian ZM, et al. (2004) Soudi China Sea 
surface water evolution over the last 12 Myr: A south-north comparison from 
Ocean DriUing Program Sites 1143 and 1146. Paleoceanography 19: PA1009, 
doi: 10. 1029/2003PA000906. 

107. Raymo ME, Rudchman WF, Froelieh PN (1988) Influence of late Cenozoic 
mountain building on ocean geochemical cycles. Geology 16: 649—653. 

108. Raymo ME, Ruddiman WF (1992) Tectonic forcing of late Cenozoic climate. 
Nature 359: 117-122. 



PLCS ONE I www.plosone.org 



12 



April 2014 | Volume 9 | Issue 4 | e95754 



